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Abstract 

Parametric energy-level correlation describes the response of the 
energy-level statistics to an external parameter such as the magnetic 
field. Using semiclassical periodic-orbit theory for a chaotic system, 
we evaluate the parametric energy-level correlation depending on the 
magnetic field difference. The small-time expansion of the spectral 
form factor K (t) is shown to be in agreement with the prediction of 
parameter dependent random-matrix theory to all orders in r. 
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1 Introduction 



More than two decades have passed since the universal energy level statistics 
was conjectured for classically chaotic systems jT]. Spectral correlations were 
found to coincide with the predictions of Random Matrix Theory (RMT). 
If the system is time-reversal invariant, the energy- level correlation in the 
semiclassical limit is asymptotically in agreement with the eigenvalue corre- 
lation of the Gaussian Orthogonal Ensemble (GOE) of random matrices. In 
a magnetic field the time-reversal invariance is broken and the energy-level 
statistics is then qualitatively affected. In that case the Gaussian Unitary 
Ensemble (GUE) gives a precise prediction for the asymptotic behavior of 
the energy-level correlation. 

Much effort has been paid to explain the agreement with RMT in terms 
of the semiclassical periodic-orbit theory|2 r A typical physical quantity, 
the spectral form factor K(t), can be written as a sum over periodic-orbit 
pairs. Berry calculated the leading contribution, of first order in the time 
variable r, by means of the diagonal approximation [Hj which is applied to 
both of the GOE and GUE universality classes. For a system with time- 
reversal invariance, the pairs of identical orbits and the pairs of mutually 
time reversed orbits both contribute to the first-order term. For a system 
without time-reversal invariance, we need to care only about the pairs of 
identical orbits. In this way one is able to partially reproduce the RMT 
prediction using periodic-orbit theory. 

Berry's work was extended to the second-order term by Sieber and Richter 
(SR) who specified the family of contributing orbit pairs 4 . The possibility 
to include more complicated orbit pairs by a combinatorial method was soon 
noticed. Heusler et al. developed the analysis to the third-order termJH] and 
Miiller et al. obtained the expansion in agreement with the RMT result to 
all orders 011]. 

On the other hand, it is also conjectured that parameter-dependent ran- 
dom matrices describe the transition of level statistics within and in be- 
tween the universality classes CHI- Saito and Nagaoinj applied semiclassi- 
cal periodic-orbit theory to the parametric transition between the GOE and 
GUE universality classes and obtained agreement with "parametric" RMT 
up to the third order. In this paper, we deal with the parametric transition 
within the GUE symmetry class, employing the magnetic field as the param- 
eter. Using semiclassical periodic-orbit theory, we evaluate the small-time 
expansion of the spectral form factor for the parametric correlation. The 
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agreement with parametric RMT is established to all orders. 

This paper is organized as follows. In § 2, a parametric random-matrix 
theory is developed and an RMT prediction for the spectral form factor is 
deduced. In § 3 and § 4, we employ periodic-orbit theory for a chaotic system 
in a magnetic field to show that a small-time expansion of the form factor 
agrees with the RMT prediction. In § 5, the key identity (a sum formula) 
used in § 4 is proved. In addition, a similar description of the GOE to GUE 
transition is briefly given in last section. 



2 Parametric Random Matrix Theory 

A parameter-dependent random-matrix theory (matrix Brownian-motion model) 
was first formulated by Dyson He considered an ensemble of iV x N her- 
mitian random matrices H which are close to an "unperturbed" hermitian 
matrix H^°\ The conditional probability distribution function of H is given 
by 

Tr{(H -e^H^) 2 }' 



P(H;a\H w ) dH (x exp 



~2a 



dH 



with 



N N 

dH = IT dH n IT dReify dhnHji. 



(2.1) 



(2.2) 



The parametric motion of the matrix H depending on the fictitious time 
parameter a is of interest. At the initial time a = 0, H is equated with 
the hermitian matrix H^°\ In the limit o — > oo, the probability distribution 
function (p.d.f.) of H becomes that of the GUE 



P(H; oo\HW) dH oc e~ TrH dH, 



(2.3) 



which is independent of H^°\ 

Let us denote the eigenvalues of the hermitian matrices H and H^ as 

,x$, respectively. Then the p.d.f. of the 



Xi,X 2 , 



, x and x^ , x^ , 



eigenvalues of if at a (under the condition that Xj 
at a — 0) can be derived as 



x 
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« n e-torw?*?'* n -13) de %fe> 4 0) )]i,^ > 2,- > iv n 

(2.4) 



where 

^ y ) = e -^ 2 ^ 2 )/ 2 £ Hj (j}^ y h -V+(V*))° (2.5) 



#,(x) = (-l)V ' . (2.6) 



with the Hermite polynomials 

_cP_ 

In the limit o — > oo, this p.d.f. becomes the p.d.f. of the GUE eigenvalues 

as 

p(xi, x 2 , ■ ■ ■ , xat; oo|4 0) , %2, ■ ■ ■ , x {0) ) = Pgve(xi,x 2 , x N ), (2.7) 
where 

N N 

Pgve(xi,x 2 ,- --jXjv) oc JJe~ (aii) ni x i _x '| 2 > ( 2 - 8 ) 

as expected. 

Now we suppose that the initial matrix is a GUE random matrix, 
so that the p.d.f. of x\ , x 2 , • ■ • , x^ is also given by (|2.8j) . Then the tran- 
sition within the GUE symmetry class (the GUE to GUE transition) is ob- 
served. The dynamical (density-density) correlation function which describes 
the GUE to GUE transition is defined as 

p d (x;a\y)=N* I -^±, (2.9) 
-<o 

where 

I(xi, o r |xj°' ) ) 

oo roo roo roc , Q . poo , q . roo . Q . 

dx? / dx* • ■ • / dxjv dxn / dx, ■ • ■ / dx, 



OO J — oo 



;at / ux 2 / ux 3 • • • / aX N 

J — oo J —oo J —oo 

I I (0) (0) (OK / (0) (0) (0)\ 

x p(xt, x 2 , ■ ■ ■ , xtv; cr|xi ; , x 2 V ■ ■ , xV)Pgue(^i , a?2 , " " " > X V) 

(2.10) 

and _ _ 

/OO /'OO 
dx/ dy/(z;a|y). (2.11) 
-oo 7 — oo 



The dynamical correlation function describes correlations between the spec- 
tra of H and H . 

It is possible to evaluate the asymptotic limit N —>■ oo of the dynamical 
correlation function |13| I14j. Introducing scaled parameters rj, X, Y as 

a = ri/(An 2 p 2 ), x = V2Nz + (X/p) , y = V2Nz + (Y/p) (2.12) 

(p = \J2N(1 — z 2 )/tt is the asymptotic eigenvalue density at V2Nz, — 1 < 
z < 1), we find 

Pd(x;a\y) 



P 2 



1 ~ p(£- n) = I du e u2 " /4 cos(7r<) f°° dv e~ v ^ /A cos(7r<), 
Jo Ji 

(2.13) 

where £ = X — Y . The Fourier transform K rm (t) = J!^^ e t2nT i p(£;rj) is 
called the form factor. For times in the interval < r < 1 the form factor 
can be written as 

K RM (r) = \ I' d M e- A ^) = ^sinh(Ar), (2.14) 
A = t]t; (2.15) 

the variable A was introduced here because it is the expansion of K-&m(t) in 
powers of r at fixed A which is most naturally connected with the semiclas- 
sical periodic-orbit theory; this expansion 

^MM=re-g^ (2.16) 

will be compared with a semiclassical result. For that purpose, we write the 
expansion into the form 

K RM (r) = Kgf(r) + X&(t) (2.17) 

with 

^(r) = re~\ K$${t) = re"* £ (2.18) 

In § 3, we evaluate the semiclassical form factor for a chaotic system and 
obtain the first-order term in agreement with K^ s \t). Moreover, in § 4, 
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the semiclassical calculation is extended to yield a result in agreement with 
the Laplace transform (taken for fixed i], using (|2.15Jl ) 



-q\ -^RM l T J 

2 



dA 



E 



r=X/rj 



i(2j + l)\J 



-(9+l)A 



A 



2i-i 



E 



1 (4j-l)! 1 



We thus show the agreement up to all orders. 



A 2 MA 
(2.19) 



3 Periodic-Orbit Theory for a Chaotic Sys- 
tem 

We consider a bounded quantum system with / degrees of freedom in a mag- 
netic field B, assuming that the corresponding classical dynamics is chaotic. 
Let us denote the energy by E and each phase space point by a 2f dimen- 
sional vector x = (q, p), where / dimensional vectors q and p specify the 
position and momentum, respectively. In the semiclassical limit H — > 0, the 
energy-level density p(E; B) can be written in the form 

p{E-B) ~ Pav (E)+ PoBC (E;B). (3.1) 

Here p scv (E) is the local average of the level density and p osc (E; B) describes 
the fluctuation around the average. 

The local average of the level density is equal to the number of Planck 
cells inside the energy shell 

where Q(E) is the volume of the energy shell. We assume that the magnetic 
field is sufficiently weak such that the cyclotron radius is much larger than the 
system size and thus the presence of the magnetic field does not significantly 
change Q(E). 

On the other hand, the fluctuation part is given by a sum over the classical 
periodic orbits 7 as 

p osc (E;B) = -±-ReJ2A,e^ E ^ B »/\ (3.3) 
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where S 1 is the classical action and A 1 is the stability amplitude (including 
the Maslov phase). The phase 7 (B) is a function of the magnetic field and 
is defined as 

6,(B) = B J a(q) • dq = B j g 7 (t) dt, g 7 (t) = a(q 7 ) • (3.4) 

where a(q) is the gauge potential which generates the unit magnetic field 
and q 7 (i) describes a classical motion in the configuration space along the 
orbit 7. 

In analogy with ([2. 13)1 . we introduce the scaled parametric correlation 
function as 

R{s;B,B') = ( P ( E + 1 ^ yi ^ P ^ ~ ^ 



Pav(^) 2 / 
Pose (-S + 2p av S (£;) 5 5 ) P°sc (-E - 2p J(E) ' - B ' 



Pav(£) S 



(3.5) 



Here the angular bracket means two averages, one over the center energy E 
and one over a time interval much smaller than the Heisenberg time 

T H = 2nhp av (E) = j^jzi ■ (3.6) 

The form factor, namely the Fourier transform of R(s; B, B'), is then written 



as 



/CO 
ds e l2nrs R(s; B, B') 
-00 

. T .o osc (E + ^B)p osc (E-^B')\ 
de e ttTT « /h - — ^ ^7—^ '- ). 3.7 

Pav(^) / 1 ' 

Putting ()3.3|) into ()3.7jl . we find that the form factor is expressed as a double 
sum over periodic orbits 

K( r ) ~ ^ /E A 7 A;,e^-^/V^-M B '))/^ ( T - ) (3.8) 



(an asterisk means a complex conjugate), where T 7 and Ty are the periods 
of the periodic orbit 7 and its partner 7', which "feel" the magnetic fields B 



7 



and B', respectively. We assume that the difference between these fields is 
sufficiently small so that its influence on the classical motion can be neglected; 
we only have to keep the resulting difference between the magnetic phases 

Let us now denote by 73- a stretch of the periodic orbit 7 whose dura- 
tion T is much larger than all classical correlation times; this stretch can 
coincide with the whole orbit (and then T is the orbit period). For times 
large compared to the classical scales mentioned, successive changes of the 
velocity dq 7 /dt can be regarded as independent random events [To], so that 
a replacement of <7 7 (i) by Gaussian white noise is justified. An average of a 
functional F[g lr ] over Gaussian white noise is evaluated as 



(( F[9,r) » 



J Vg 1 exp 








J Vg„, exp 







(3.9) 



and implies a correlation ((g 7 (t)g 7 (t'))) = 2D5(t — t'). Including this Gaus- 
sian average (carried over the whole duration of the periodic orbits), we 
rewrite the form factor as 

H \ 7 , 7 ' H I 

(3.10) 

We shall evaluate the small-r expansion of this semiclassical form factor, 
restricting ourselves to homogeneously hyperbolic systems with two degrees 
of freedom (/ = 2). 

Let us begin with adapting Berry's diagonal approximation j3] to correla- 
tions between two spectra pertaining to different values of the magnetic field. 
In this approximation, one first considers the contributions of periodic-orbit 
pairs 7' = 7. The key ingredient is Hannay and Ozorio de Almeida (HOdA)'s 
sum rule [To] 

^£IA| 2 <^-^)=t. (3.11) 

Using this sum rule and the Gaussian average ()3.9[) for pairs of identical 
orbits (7, 7) we find 

^EIV<^ - ^) <<e^)-W)A>> = re-". (3.12) 

J-H 7 \ 1 H/ 
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Here T is the period tT#. Since the Heisenberg time T# is of the order 
1/h and a = (B — B') 2 D/h 2 the decay rate at r fixed is proportional to 
(B - B') 2 /h 3 . The contribution of pairs of identical orbits does not vanish in 
the limit h — > provided the field difference is scaled such that this parameter 
remains finite. 

Consider now the case when the magnetic field is so weak that its influence 
on the orbital motion can be neglected. Then the system is close to being 
time-reversal invariant, and its periodic orbits occur in almost mutually time- 
reversed pairs (7,7); these must be taken into account as well. However we 
can check that the pair (7, 7) yields no contribution. This will be true if both 
B and B' are quantum mechanically large in the sense 

B,B' > 0{h 3/2 ), (3.13) 

which does not prevent the field difference from being quantum mechanically 
small. Namely, as the phase factor 9 1 changes sign under time reversal, 

-\T \A,\ 2 5 (t - ((e W)-<W)/ft)) = r((e W)+W)A)) 

(3.14) 

in the limit h — > 0. It means that pairs of time reversed orbits do not 
contribute to the form factor if ()3.13j) holds. 

Putting the above results together, we obtain the diagonal approximation 
of the form factor as 

KpQ g ^ = re~ aT . (3.15) 

This is in agreement with the first-order term of the RMT prediction ()2.16|) . 
if the RMT parameter A is identified with aT. 

4 Off-diagonal Contributions 

We are now in a position to calculate the off-diagonal contribution. In order 
to identify the family of periodic-orbit pairs responsible for the leading off- 
diagonal terms, we note the fact that long periodic orbits have close self- 
encounters where two or more orbit segments come close in phase space. 
The duration of the relevant self-encounters are of the order of the Ehrenfest 
time 7 / J7| . Although T# is logarithmically divergent in the limit h — > 0, 
it is still vanishingly small compared to the period (which is of the order 
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of the Heisenberg time Tjj). After leaving a self-encounter, the orbit goes 
along a loop in phase space and comes to a different (or back to the same) 
encounter. All off-diagonal terms arise from the existence of orbits 7 which 
are close but different from the partners 7' in the encounters but almost 
identical to them on the loops. Within the encounters the orbits 7 and 7' 
are differently connected to the loops. Suppose that the magnetic fields B 
and B' are sufficiently strong. Then, since we are treating a system without 
time-reversal invariance, 7 and its partner 7' go in the same direction on all 
loops. 

Let us consider such a periodic-orbit pair a = (7,7') with L loops and 
V encounters. Inside each encounter, we introduce a Poincare section V 
transversal to the orbit 7 in phase space. Pairwise normalized vectors e s and 
e u span the section V. Here the vectors e s and e u have directions along the 
stable and unstable manifolds, respectively. Each segment of the orbit within 
the encounter pierces through V at one phase-space point. The displacement 
5x between such points can be decomposed as 5x = se s + ue u . If we fix 
one reference piercing point as the origin, each of the others is specified by a 
coordinate pair (s,u). 

Suppose that the periodic orbit 7 pierces V within the r-th encounter. If 
l r segments of 7 are contained in the encounter, there are l r piercing points 
so that l r — 1 coordinate pairs relative to the reference piercing are necessary 
to specify them. Consequently, we need J^=i(^r — 1) = L — V coordinate 
pairs (sj,Uj) to specify all the piercing points within the encounters. 

We denote the time elapsed on the j'th loop by Tj and the duration of the 
r-th encounter by t enCjr . It follows that the total duration of the encounters 
is 

v 

t a = ^ ^ lrtenc,r • (4. 1 ) 

r=l 

Using these notations, we can employ ergodicity to estimate the number of 
encounters in a periodic orbit with a period T = Y%=i Tj + t a as jHl 13 El HI] 

rT—ta rT—t a —Ti rT—t a —Ti-T2 Tl-2 

duds/ dTi / dT 2 .../ dT L ^Q a , (4.2) 



Jo JO 



where 



Qa = N(v) — T (4.3) 

llr=l (/ enc,r 1L 
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and the integration measures are given by 

L-V L-V 

du = n dn j' ds = n ds i- ( 4 - 4 ) 

i=i 3=1 
The combinatorial factor N(v) is the number of structures of orbit pairs for a 
given vector v = (v 2 , v 3 , t> 4 , • • •), where the component v\ denotes the number 
of the encounters with I segments; we will occasionally write 

v = (2) 1 ' 2 (3) U3 (4)" 4 (4.5) 

It should be noted that 



oo oo 



L = X>,, V = '£v l . (4.6) 

1=2 1=2 

For n — L — V + l — 3 and 5, we tabulate N(v)'8 in Table 1. The precise 
meaning of the orbit structure is expounded in next section. 

We then calculate the Gaussian average (|3.9|) on the loops and obtain 
a factor e _aTl e _aT2 ■ • • e~ aTh . Similarly, an encounter contributes a factor 

It is now straightforward to obtain the contribution to the form factor 
from the orbit pair a 



KpoA T ) = T J duds 



x 



T— t a pT—t a —Ti /■T—t a -T 1 -T 2 T L _ 2 . A rf/« 

m dT 2 --- dT L ^Q a R a e tAS /\ 



o Jo Jo 



(4.7) 

where 

fi^ = e -a(T 1 +T2+---+T L ) e -a((l 1 ) 2 t cnc , 1 + (l 2 ) 2 t cnc ,2+---+(lv) 2 t C nc,v) U g\ 

The action difference AS = S^ — Sy is estimated as AS = J2j=i u j s j [HHZ1IH1 
This formula contributes to the terms of order r n with n = L — V + 1. 

Then we expand Kp 0ta (r) in t enCjr and extract the term where all t e nc,r's 
mutually cancel. Because of the appearances of extra factors h or rapid 
oscillations in the limit h — > 0, the other terms give no contribution [HI IE]- 
We thus obtain the off-diagonal term of the form factor 

<V) = E^x(^)^ V_1 fl(-^|;-W 2 a)/(T) 

- X>w T y n (-'„-''«) nn (4.9) 
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where 

/•T (T-T\ rT-T 1 -T 2 T L _ 



f{T)= [ dTi / 1 dT 2 ---( 1 2 i_2 dT L _ 1 e- a(Tl+T2+ - Ti) . (4.10) 
Jo Jo Jo 



Let us put A = aT and calculate the Laplace transform of K^q\t)/t 2 as 
f e-^^-^dA 

Jo T 2 

1 / 1 \ L ~ V - 1 roo 00 / 8 \ Vl 

_ / i s L—V—l oo i 



oo i / 1 \ n ~ 2 L-V+l=n oo / i \ 

(4.11) 

where N(v) = N(v)(—l) v H1Z2 l Vl /L. In the above equation, a simple graph- 
ical rule is observed: each loop contributes a factor l/(a(q + 1)) and each 
encounter contributes —la(q + I). In next section, we shall prove a sum 
formula for n > 2 

n-l 



L-V+l=n 00 / 1 \ 



(2n - 3)! / 1 



1 \ - , n odd, 



(J, n even, 

(4.12) 

from which it follows that 

Jo t' ,t;(or H )^i(2j + l)!(9 + l)«' 1 J 

As aTff = (aT)(Tff/T) = X/r = rj, this is in agreement with the RMT result 

Ann. 



5 A Sum Formula for N(v 



In this section we shall give a proof for the sum formula (see (|4.12j) ) 
L-v+i=n 00 f (2n — 3)! n _ x 

E #00 n(l + (/ - = \ n\ X ' n ° dd ' (5.1) 
« 1=2 I 0, n even 
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with n > 2. For that purpose we introduce a number N P (v) depending on 
the vector 

# = ^yi (2)^(3)^3 (4)«4... (5.2) 

and set L = ^ an d V = X^i Let us denote an "encounter" permu- 
tation of the numbers 1, 2, • • • , L as 

PCI1C = ( Pc„c(l) Pcnc(2) P cn !(3) ••• P^L)) ^ 

and define a "loop" permutation 

( 1 2 3 ••• L — l L \ 
P ^={ 2 3 4 -.. L 1 J' ^ 

We define Np(v) as the number of permutations P cnc which satisfy the fol- 
lowing two conditions. 

(A) The permutation P enc has vi cycles of length I. 

(B) The product Pi 00 pPenc is a permutation with a single cycle. 
Then it follows that 

N((2) V2 (3) V3 (A) V4 ■■■) = iVp((l) (2) ^,2 (3) ^,;i (4) ^,4 • • •). (5.5) 

In order to explain the reason, let us suppose the following situation. The 
encounters include Y^=\ l r = L orbit segments in total, so that there are L 
"entrances" where the orbits come in and L "exits" where the orbits go out. 
A periodic orbit 7 comes in an encounter at the first "entrance" and goes 
out at the first "exit" . Then it comes to the second "entrance" and goes out 
at the second "exit" . It continues to follow the connection pattern 

j-th "entrance" — > j-th "exit" — > (j + l)-th "entrance" 

and finally goes out at the L-th "exit" and then comes back to the first 
"entrance" again. On the other hand, the partner orbit 7' comes in an 
encounter at the first "entrance" and goes out at P cnc (l)-th "exit". Then it 
must go to the Pi 00 pPenc(l)-th "entrance", as the partners go along the same 
loop. It continues to follow the pattern 

j-th "entrance" — > P enc (j)-th "exit" — * Pi 00 pPenc(j)-th "entrance" 

In this manner, if a permutation P enc is given, the structure of a periodic 
orbit 7' is specified. 
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The j-th "entrance" and the l-ih "exit" belong to the same encounter, 
if and only if j and I are contained in the same cycle of the permutation 
Penc- Hence the condition (A) is required. The orbit 7' finally comes to 
(Pi oop P enc ) L (l)-th "entrance". As 7' is a connected periodic orbit, it must be 
the first return to the first "entrance" . This is guaranteed by the condition 
(B). 

A combinatorial argument [HI El IE] yields a recursion relation for 

00 

N P (v) =N P (v)(-l) v l[ l*/L (5.6) 

1=1 

as 

Vl N P (v) + Y,v^rT l ' 1] kN P (^ k+l - 1] ) 

k>l 

+ E {vi- m -x + l)vt* m ' l ~ m ~ l] Np{^ m ' l ~ m ~ l] ) = ^ (5.7) 

\<m<l-2 

Here we used a notation 

^ l a.i,-,a.v->l3i,-,P v i] ^ (5-8) 

which is the vector obtained from v when we decrease each of v ai , v a2 , • • • , v av 
by one and increase each of vp 1 , Vp 2 , ■ • • , Vp u , by one. It should be noted that 
Np{v) is zero if any of the components of v is negative. 

In the special case / = 2, we obtain a simplified recursion formula for 

N{v) 

v 2 N(v) + £ v [ ^ k+1] kN(v^ k+ V) = 0. (5.9) 

k>2 

Let us introduce a variable x and define 

00 

N(v,x) = N(v)H(l + {I -l)x) Vl . (5.10) 
1=2 

Then the recursion formula (J5.9)) reads 

») + E fetl t ^ ~ 1)a:) 4yr fc+ '^(^' w ". *> = o- (5.H) 

Summing this over u with fixed L — V + 1 = n, we find 



L-V+l=n 

E 



\ -\- X k>2 ^ 



= 0. 
(5.12) 
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Here the sum over v can be replaced by the sum over v' = ■y[ fc > 2 ^ fc + 1 l ) so that 

L-V+l=n 

E i^i^s). (5.13) 



L-V+l=n 



E * 2 r fc+1] iv(#' 2 - fc+i U) 



Dropping the primes, we can thus write 



L-V+l=n 

E 

if 

L-y+l=n 

E 



1 + x 



w 2 „ fc(l + (k - l)x) 

+ E Y+fac — Vk+l 

k>2 

J2 v k (k - 1) - E - 1 " >x 

k>2 k>2 1 + - 1)^ 



7V(tT, x) 



n — 1 — x 



dx 



L-V+l=n 

E iV(i?, ar) = 0, 



(5.14) 



which means 



L-V+l=n L-y+l=n 

E N(v,x)=C n x n -\ c n = e 



(5.15) 



Thus the sum formula has been proved up to a constant C n . 

Let us then calculate C n . First note that, according to (|5.10p . each N(v, x) 
contains only terms of the order x v and lower orders. Due to the inequality 



n - 1 - V = L - 2V = E M l ~ 2) > 0, 

1=2 



(5.16) 



this means that the largest order possible for a given n = L — V + lis x n_1 . 
This order is reached only for v with v% = v± = ■ ■ ■ = 0, for which the equality 
holds in (|5.16j) . Accordingly, we find 



L-V+l=n 

E N(v,x) 



L-V+l=n 



E N(v)x V Il(l-l + -) 

v 1=2 V X/ 

= A^((2) n ~ 1 )x n_1 + lower order terms in x. (5.17) 
Comparison with (|5.14|) now yields 

C n = N((2) n - 1 ); (5.18) 
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all terms of lower orders in x must mutually cancel. In order to evaluate 
A^((2)™ _1 ), we can utilize a closed expression for Np(v) (with Vj > for 
j < A and v j = for j > A) 

A A 

= EE-E (-D^-»^ J =' A g , 



nM«i-^)q 

J'=l 

(5.19) 



which was derived by Jiirgen Muller^Tj. Using the identity 

e - s s j ds = j!, (5.20) 
we can rewrite Jiirgen Miiller's formula as 

f—l)V roo roo oo (njj _ (_ x )j)vj 

Np(v) = -^- dx dye-e-*!! - , . ( 5 - 21 ) 
L[L + 1) Jo Jo 7_i Vi\ 



3=1 "3- 



so that 

A>((2) n - 1 ) = A> P ((2) n - 1 ) 



( 1 ^tt— 1 />oo /"oo 

/ dx / dy e- x e- y {y 2 -x 2 ) 71 - 1 
Jo Jo 

ds f dt e-'s n -H n - 1 



2(n-l)(2ri-l)(n-l)!io ' Vo 
(-I)"' 1 



4(n-l)(2n-l)(n-l)!io 
(2n-3)! 

n! ' nodd ' (5.22) 
0, n even 

(s = x + y,t = x — y), which establishes the desired result (|5.ip . 

It is easy to check that Jiirgen Miiller's formula holds for if s with small 
L — V (for example, A r p((l) 1 ) = —1). Therefore, in order to prove it in 
general, it is sufficient to verify that it fulfills the recursion relation ([5.7)1 . 
For that purpose, we first define an "average" (■ • -) tt of a function f(x, y) as 

roo roo °° 

(f(x,y)),= dx dye- x e-vf(x,y)l[(yi -(-xyp. (5.23) 
Jo Jo t=i 
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Since Np(v) = if any of Vj is negative, (|5.7|) evidently holds if v i = 0. Hence 
we focus on the case vi > 1. Then partial integrations yield a relation 



y l 1 + (— x) 
y l - (-x) 1 



i-i 



d 



d 



,21-1 



<!>*-* 

^* ~{dy \y l - {-x) 1 J dx \y l - {-x) 1 J ' (y l - (-x) i y 

roc roo 

/ dx / dye~ x e~ y 
Jo Jo 

y i q ^ (_ x y 2i-i _ (_ x )H-i 
y l — {—xf dy y l — (—x) 1 dx (y l — (— x) 1 ) 2 



\2l-l 



3=1 



(5.24) 



for I = 1, 2, • • • , L. Using this relation and the identity 

y l d (—x) 1 d 



y l — (—x) 1 dy y l — (—x) 1 dx 



J 3=1 



-X 



yk+i-i _ r_ x \k+i-i x 



lk kVk '(v l - (-*)W - (-*)*) fJi 



(5.25) 



we can readily derive 



y 



i-i 



-X 



,1-1 



y l - (-x) 1 



hi) 



k>l 



The following identity 



y k+i-i _ (_ z )fc+i-i \ 
{y i _ (- X )i)( y k _ / 9 - 

(5.26) 



OO POO 



1 OO 

dx I dye™ e-^—-Jl(yi 
x + y 



JO 
-L-l 



roo roo I °° 

/ dx dy e- x e- y —- J[(y j 
Jo Jo x + y =1 



\3\ v i 



-x 



\3\ v 3 



(5.27) 



can be proved by a transformation of the variables ujx \— » x, uy \— > y. Differ- 
entiating the both sides of this identity with respect to u and then putting 
u — 1, we obtain a relation 



1 



Y <i>* 



l ' 



(5.28) 
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from which it follows that 



L + l x ,v \y l -(-x) 1 [" x + y 

Then, utilizing 



xy l 1 + y(—x) 1 1 1 



- £ [("*) "V - ^ 1 + {-x) l - m - x y m ] , (5.30) 



we find 



2 /^i + 



L + l Wt ' \ y'-(-x)' I, 

/ ( y ™ - (-x)" 1 )^-"- 1 - (-x) l - m - 1 ) \ (53r 
i<m<z-2 \ y l ~ (~ X Y I $ 

Adding the both sides of ()5.26|) and f)5.31j) . we arrive at 

T 1 / „fc+J-l /■ -U+J-l \ 

^ <!>*=£* ((«*-*«) 
^ + 1 fc>l \ 



+ £ / ( f -(-^)(^-i-(-^-i) \ ^ 

l<m<Z-2 \ V 1 ~ ( — x )' / 

which gives the desired recursion relation ()5.7|) with Jiirgen Miiller's formula 
dnHH) substituted. 



6 The GOE to GUE Transition 

The equal-parameter correlation function i?(s; B, B) describes the transition 
between the GOE and GUE universality classes as the magnetic field B 
increases from zero [TTJ HBJ CHI- I n this section, we shall reproduce Saito and 
Nagao's semiclassical calculation^!] of the form factor (the Fourier transform 
of R(s;B,B)) and further derive a sum formula analogous to ()5.1|) as a 
conjecture. 

The RMT prediction of the form factor in this case is derived from Pandey 
and Mehta's two-matrix model[20j. For small r (0 < r < 1), it can be written 
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as 



2 ii-2r k + 2r 



o 1 + 2ry 

(6.1) 

In the GOE limit the parameter fj, is zero and in the GUE limit it goes to 
infinity. 

The semiclassical argument is similar to that in § 3 and § 4. The difference 
is that we have to take account of the mutually time reversed pairs of loops 
and segments of classical orbits. Following a similar argument as in § 3, we 
obtain a diagonal approximation for the form factor 

4o aS) (r) = r + re~ bT , (6.2) 
b = 4B 2 D/h 2 . (6.3) 

The RMT parameter \x should be equated with bT in reference to the semi- 
classical result. 

In order to extend the calculation to the off-diagonal terms, we need 
to introduce integers n enCir and M characterizing the structure of the orbit 
pairs as follows. Let us fix an arbitrary direction (+) in which the orbits pass 
through the r-th encounter and call the opposite direction (— ). Suppose that 
the orbit 7 passes through the encounter #^(7) and #( - )(7) times in (+) 
and (— ) directions, respectively. We then define the number n en c,r as 



1 



{# (+) (7) " # H (7)} " {# (+) (7) " # ( " ) (7 / )}| • (6-4) 



Moreover we define M as the number of the pairs of mutually time reversed 
loops. 

As before, for a general orbit pair a with L loops and V encounters, the 
number of encounters in one periodic orbit of a period T is evaluated as 

/rT—ta. rT—ta-Ti rT—t a -T\—T2 Tl-i 

duds J o dT x J q dT 2 --- ^ dT L _ x Q a , (6.5) 

where 

Q a = N(v,M) T (6.6) 

^ 1 lr=l "enc,r J ' 
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Here the combinatorial factor N(v, M) depends on a matrix v and M. The 
component V\ m of the matrix v is the number of the encounters with l r = I 
and n enc>r = m. One can write 

v = (2, 0) V20 (2, If 21 (2, 2) V22 (6.7) 

Following the argument in 00|B], we can identify N(v, M) with the number 
of generalized permutations satisfying suitable conditions. 

Let us consider the effect of the gauge potential. The Gaussian average 
(|3.9jl on the loops gives a factor e _bTl e~ 6T2 • ■ ■ e~ bTM , while from an encounter 
it yields e~ 6 ^ ncnc - r ^ * cnc r . Thus we conclude that the total contribution to the 
form factor from the orbit pair a is 



Kpo,a( T ) — t J duds 



X 



T— t a pT—t a —Ti rT—t a —T\—T2 Tl-i . «r./» 

dT x \ dT 2 --- dT L ^Q Q R Q e lAS/h 



o Jo Jo 



(6.8) 

with 

= e ~K T \+ T 2-\ l-r M ) e -6((n cnCi i) 2 t cnc ,i+(?i cnCi 2) 2 tonc,2H h(n cnCiV -) 2 t C nc,v) 

This contributes to the terms of order r n with n = L — V + 1. As before 
we expand Kp 0tQ (r) in t enCtr and extract the term where all t e nc,r's mutually 
cancel. Then we find that the off diagonal contribution to the form factor is 



V M=0 L KI HS l=2m=0\ 01 J 

(6.10) 

where 



r-T r-T — Ti rT — Ti — T2 Tl — 2 

f(T,M)= dTi / dT 2 ---/ dT L _ ie - b(Tl+T2+ - TM) . 



Jo Jo 

(6.11) 

If M = 0, the direction of motion along all loops and hence in all encounters 
does not change in the partner orbit; consequently n enC)T = for all encoun- 
ters. The corresponding structures also exist in the case without time-reversal 
invariance, so that 

N(v o) = J N((2)^(3)^ ■ ■ •), if all v nj with j ? vanish, 
|0, otherwise. 
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Here N(v) is the number of structures introduced in § 4 and § 5. Time 
reversal of each such partner orbit produces another partner with M = L; 
therefore 

N(v L) = { N (( 2 ) V22 (3) V33 ■ ■ •)> if all v nj with j ± n vanish, 
1 0, otherwise. 

Note that the structures with M = 0, L may exist only for odd n = L — V+l; 
see 0. 

Noting the above relations for the combinatorial factors, we can eval- 
uate the contribution of the structures with M = 0, L in the same way 
as in § 4, namely by Laplace-transforming the corresponding summands in 
Kpq\t)/t 2 , using the sum rule (j5.1j) for N(v) and transforming back to the 
time representation. The contribution of the structures with M = turns out 
to be zero whereas the structures with M = L reproduce the third summand 
in the last line of (JbMj) . On the other hand, making the Laplace transform of 
the part of KpQ \t)/t 2 with 1 < M < L — 1 and equating the result to the 
corresponding RMT prediction deduced from the integral in the last line of 
()6.ip . we arrive at a conjecture 

L-V+l=n L-l N S M S ( X y oo oo 
V M=l ^ l iT ^ l=2m=0 

1 n_1 / x \ n -p- 1 n_1 (—l) j 2 j 

(l + ^SvTT^J (2n-p-3)!g i(n _ i _ 1)!(j _^ )r 

(6.14) 

In the cases n = 2 and 3, the conjecture (|6.14j) was substantially proved in 
[TTj ; by machine-assisted counting it was verified up to n = 7. For small val- 
ues of n up to 4, the relevant N(v, M)'s are tabulated in Table 2. Moreover, 
putting x = 0, we obtain 

TET'-'i'n^i-r^ (6.i5, 

v a/=i ■ L 1=2 n i 

(vi = J2m=o v im), which is relevant to the GOE form factor. This special 
case was proved in |S1 13 El- The full proof of (jbM4j) is an interesting open 
problem. 
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n 


V 


L 


V 


iV(V) 


3 




4 


2 


1 






3 


1 


1 


5 


(2) 4 


8 


4 


21 






7 


3 


49 






6 


2 


24 




(3) a 


6 


2 


12 




(5) 1 


5 


1 


8 



Table 1: The number N(v) of the orbit structures corresponding to the vector 
v= (2) V ^(3) V3 (A) V4 L = Eilvi, V = and n = L -V + I. 
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n 


V 


L 


V 


M 


jV(v, M) 


2 




2 


1 


1 


2 


3 


(2,0) 2 


4 


2 


2 


4 




(2,0) 2 


4 


2 





1 




(2, 1) 2 


4 


2 


2 


4 




(2,2)* 


4 


2 


4 


1 




(3, or 


3 


1 


2 


3 




(3, or 


3 


1 





1 




(3, I) 1 


3 


1 


1 


3 




(3,sr 


3 


1 


3 


1 


4 


(2,0) 3 


6 


3 


1 


6 




(2, or 


6 


3 


3 


10 




(2, 0r(2, l) 1 


6 


3 


2 


12 




(2, or (2, ir 


6 


3 


3 


36 




(2,0) 1 (2,l) i (2,2) 1 


6 


3 


4 


12 




(2,2) a (2,0) 1 


6 


3 


5 


6 




(2,0) 1 (3,0) 1 


5 


2 


1 


15 




(2, (J) 1 (3, 0) 1 


5 


2 


3 


15 




(2, 0r(3, l) 1 


5 


2 


2 


25 




(2,0) i (3,2) 1 


5 


2 


3 


10 




(2,0) 1 (3,3) 1 


5 


2 


4 


5 




(2,l) i (3,0) i 


5 


2 


2 


20 




(2,l) 1 (3,l) i 


5 


2 


3 


20 




(2,2) A (3,1) A 


5 


2 


4 


10 




(4,0)' 


4 


1 


1 


12 




(4,0)' 


4 


1 


3 


4 




(4, 1) 1 


4 


1 


2 


16 




(4,2)' 


4 


1 


3 


8 



Table 2: The number N(v,M) of the orbit structures corresponding to the 
matrix v = (2, 0)^(2, 1)^(2, 2)™ •• •, L = £i£r>Jm, V = J2iJ2 m vi m , 
n = L — V + 1 and the number M of the pairs of mutually time reversed 
loops. By machine-assisted counting the table is extended to higher values 
of n. 
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